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ABSTRACT 


An experimental program to obtain uniform inlet flow to 
the test blading in a large cascade facility designed to use 
inlet turning vanes, and to measure the conventional blade 
element performance, is described. Attempts to reduce non- 
uniformities (+1% in velocity) using screens were unsuccess- 
ful and so abandoned. Preliminary DCA biade element 
performance data were obtained without screens at one in- 
cidence angle before aero-mechanical problems with the inlet 
guide vane assembly curtailed testing. The blade surface 
pressure distribution at the one test condition compared 
very favorably with the distribution predicted using the | 
NASA computer code QSONIC. Recommendations were made that 


wculd avoid the aero-mechanical problems encountered. 
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I, INTRODUCTION 


The need for lightweight, fuel efficient gas turbines 
that are capable of developing large amounts of thrust or 
power has motivated a continuing drive to obtain more ac- 
curate predictions of the flow through turbomachinery. 
Cascade testing of blade rows has, in the past, been a logi- 
cal and relatively inexpensive way to learn more about the 
phenomena involved in the flow through compressor and tur- 
bine stages. It is required more today in order to verify 
two dimensional and near-two dimensional analysis codes for 
flow through cascades. Such testing also provides two- 

we dimensional blade element performance data which, in the 
absence of reliable analytical predictions, are required in 
the design of compressors and turbine stages. Reference 1 
describes how cascade measurements are obtained using a cas- 
cade wind tunnel and then used in the design process. 


Before subsonic cascade wind tunnel data can be accepted 


as being valid, the flow conditions must meet three require- 
ments. These criteria are discussed in detail in Refs. l, 2 
and 3. First, any disturbance in the airflow should be 
‘es caused by the test blades; that is, the inlet flow to the 
test section must be acceptably uniform. 

Secondly, the measured flow characteristics should, 


by ideally, be independent of spanwise position along the test 


15 


blades. The flow, ideally, should be two dimensional. 

Duval [Ref. 3] demonstrated that excellent flow conditions 
could be achieved in the Naval Postgraduate School Turbo- 
propulsion Laboratory (NPS/TPL) Subsonic Cascade Wind Tunnel 
using test blades with an aspect ratio of approximately two. 
The absence of suction along the walls results however, in 
some degree of streamline contraction which is measured in 
terms of an Axial Velocity Density Ratio (AVDR). 

The third requirement which must be satisfied is the 
periodicity of the inlet flow to the test section and of the 
outlet flow. Within one chord length of the leading edges 
of the test blades an upstream perturbation occurs as the 
streamlines adjust to negotiate the blade passages. Since 
the rectilinear cascade is simulating an infinite cascade of 
blades, the flow characteristics should be the same at cor- 
responding axial and blade-to-blade positions within each 
blade passage. This same condition should be true at any 


measurement plane downstream of the test blading. 


As described by Rose and Guttormsen [Ref. 4] several 
_ unique features were incorporated into the design of the 
NPS/TPL Cascade to ensure a two-dimensional and periodic 
flow at the test blading. Initial evaluations of the fa- 
e cility were conducted and reported in Refs. 3, 4, 5 and 6. 
Work by Moebius [Ref. 6] involved modifications to the tun- 


nel plenum chamber which established satisfactory uniform 
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flow at the exit of the bellmouth contraction into the test 
section. 

In order to maintain an aspect ratio close to 2.0 ata 
solidity of 1.67 Cina [Ref. 7], following the work of Duval 
(Ref. 3], used a cascade configuration of 20 blades with 3 
inch spacing. Cina conducted a program of tests of DCA 
blading at five (5) different air incidence angles. With 
this cascade configuration, Cina found that the inlet flow 
to the test section was uniform in direction and of uniform 
static pressure, but with an imposed variation in velocity 
and stagnation pressure resulting from the wakes of inlet 
guide vanes. Although excellent periodicity was found cver 
pairs of test blades, departure from strictly periodic con- 
ditions were detected from one blade passage to another. 
Cina explained this condition as being the result of the 
inlet guide vane wakes being separated at two inch intervals 
and entering a test section configured with a three inch 
blade spacing. Because of these flow conditions, Cina con- 
sidered his results to be preliminary. 

As a result of these findings the Cascade Wind Tunnel 
was modified so that inlet guide vanes were provided at one 
inch intervals. The object of the study reported herein 
was to obtain blade performance data on Cina's cascade with 
good periodic and uniform flow conditions. A necessary 
condition was to obtain agreement in the results for blade 


forces evaluated from surface pressures and from a momentum 
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balance. A second objective was to compare measured blade 
surface Mach numbers with the results of code calculations. 

At the outset, it was first necessary to carry out an 
extensive testing program to verify the new inlet guide vane 
section and the effect the new spacing had on flow uni- 
formity and periodicity. It was found that the uniformity 
of dynamic pressure was improved with the inlet guide vanes 
spaced at one inch intervals. Attempts were made to further 
improve the flow by the use of (various) wire screens placed 
downstream of the inlet guide vanes. These methods proved 
unsuccessful for the range of parameters tested and in fact 
aggravated the situation. 

Cina's testing of the Double Circular Arc blading was 
repeated without screens and with the Cascade Wind Tunnel 
configured with the modified inlet guide vane arrangement. 
Limited measurements were obtained before aero-mechanical 
problems with the new IGV arrangement, at the higher tunnel 
speeds, were encountered. 

The overall purpose of the testing program initiated by 
Cina was to obtain data with which to verify design optimi- 
zation computer codes developed by NASA. Towards this goal 
a fast, reliable computer analysis code (QSONIC) for cal- 
culating the flow field about a cascade of arbitrary 2=-D | 


airfoils was obtained from NASA. The code was adapted and 


modified to run on the Naval Postgraduate School's IBM 


370/3033 computer. 


we 


The program QSONIC was developed by NASA to overcome the 


Mach number limitations of the earlier program TSONIC (Ref. 
8].  QSONIC is described in Ref. 9. Procedures for using 
the program QSONIC at the Naval Postgraduate School's com- 
puter facility are given in Appendix D. The procedures are 
documented for the case of the DCA blading in the NPS/TPL 
cascade wind tunnel. A program listing is included to docu- 
ment changes made to the code in order to adapt to the oper- 
ating system of the NPS computer. 

Preliminary results show that experimental measurements 


and code predictions are in very good agreement. 
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II. FACILITY DESCRIPTION AND MEASUREMENT APPROACH 


A. SUBSONIC CASCADE WIND TUNNEL 

The Naval Postgraduate School's Rectilinear Cascade 
Facility is shown in Fig. 1. A description of the facility 
as it was originally configured is given in Ref. 4. The 
test facility is an open cycle wind tunnel, designed for 
the purpose of testing cascades of axial-flow turbomachinery 
compressor or turbine blades. The unique design of the test 
section ensures that the airflow paths from the inlet guide 
vanes to all of the blades of the cascade test section are 
of equal length. This particular design was intended to 
eliminate the problems found in other cascade wind tunnels 
caused by having wall boundary layers o€ different thick- 
nesses entering the cascade at different points. 

As a result of the work reported in Ref. 1, two fine 
mesh screens were installed at the bellmouth entrance to 
improve flow stability. A follow-on study into the cascade 
performance was conducted by Bartocci and is reported in 
Ref. 5. As a result of Bartocci's findings, plenum turning 
vanes were installed to direct plenum inlet air towards the 
bellmouth entrance and to decrease the total pressure fluc- 
tuations present at the bellmouth entrance. Figure 2 shows 
the configuration of the plenum chamber as modified by S8ar- 


tocci. Reference 6 describes work by Moebius that resulted 
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in further modification to the plenum chamber in which the 


original contraction was changed to two two-dimensional con- 
tractions in series. After this modification, acceptably 
small variations in velocity and flow angle were measured 
at the inlet guide vane station. Figure 3 shows the in- 


ternal arrangement of the plenum chamber as modified by 


. Moebius and as it was configured for the work presented 


here. 

. Using the plenum configuration shown in Fig. 3, Duval 
(Res 4] found that the wakes from the inlet guide vanes 
were not mixed out at the lower measuring plane of the test 
section but gave a well defined periodic variation in the 
impact pressure. The peak-to-peak variation was +4% of 
dynamic pressure over two-inch periodic intervals. This 
condition was undesirable, but was tolerated while looking 
only to establish the values of parameters required to 
achieve two-dimensionality and periodicity. Since the 
inlet flow conditions were not uniform, mass averages were 
used to calculate properties at the inlet plane from probe 
measurements. 

In order to achieve a solidity of 1.67 ane aspect ratio 
of about 2, a blade spacing of 3 inches was required for the 
tesis carried out by Cina [Ref. 7]. The tests showed unac- 
ceptable departures from blade-to-biade veriodicity under 
conditions of high blade loading and the installation of 


additional guide vanes was recommended. The modification 
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to the inlet guide vane section of the tunnel resulting from 
Cina‘'s findings is described in detail in Appendix A. 

In the present work, several tests were completed with 
the tunnel further modified by the introduction of wire 
screens between the inlet guide vanes and the lower plane of 
the test section. Appendix B describes the screen material 
and the criteria used to select the particular screens used 


in this study. 


B. INSTRUMENTATION 

The instrumentation used in the present study is that 
whica is described in detail in Ref. 7. Twenty static pres- 
sure taps were located on the north and south side walls. 
The taps on the south wall were connected to a water mano- 
meter board so that the uniformity of the static pressure 
distribution of the inlet and outlet could be monitored 
visually. Additionally, one upstream tap on each wall and 
one downstream tap on each wall (near the centerline) were 
also connected to the Scanivalve so that these static pres- 
sures were recorded. 

Figure 4 shows the probe that was used for the upstream 
survey (at the lower plane). The probe was a United Sensor 
Corporation DA 125 probe, serial number A847-1, calibrated 
earlier at various Mach numbers and yaw angles in a cali- 
bration facility. The United Sensor Corporation DC-125-24- 


F-22-CD probe, serial number A981-2 (Fig. 5), which was used 
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at the upper plane was similarly calibrated. The charac- 
teristics of the probes were approximated analytically to 
facilitate automatic data reduction procedures. The cali- 
bration and application procedures were those given by Duval 
{Ref. 3). Appendix B of Ref. 3 describes both the upstream 
and downstream probes in detail. The mounting and traversing 


mechanisms are described in Ref. 7. 


C. REFERENCE MEASUREMENTS 

Plenum chamber (supply) pressure and temperatures, and 
atmospheric pressure were recorded on each data scan. 
Plenum pressure was also displayed on a water manometer 
board. The total temperature in the test cascade was as- 


sumed to be the same as the plenum chamber temperature. 


D. TEST BLADING 

The double circular are test blading modeled the midspan 
section of the stator of the compressor stage reported in 
Ref. 10. Coordinates describing the profile of the blading 
are listed in Table D-2. The leading edge and trailing edge 
are shown in detail in Fig. 6. A photograph of the center- 
most blade is shown in Fig. 7. 

The three blades centrally located in the cascade were 
constructed with surface pressure taps along the midspan 
section as shown in Fig. 8. The centermost blade had 19 
ports on each of the pressure and suction surfaces and one 


tap at the leading edge. The two blades adjacent to the 
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center blade had 3 surface pressure taps located on each of 
the pressure and suction surfaces. The surface pressure tap 


locations for the centermost blade are given in Fig. 9. 


E. DATA ACQUISITION, REDUCTION AND ANALYSIS 

Data were recorded, reduced, and plotted using the modi- 
fied Hewlett Packard HP-3052A Data Acquisition System shown 
in Fig. 10. Reference 11 describes the system in detail. 
The system incorporated a HP-9845A desktop computer as a 
controller, with all components connected on the HP-98034A 
HP~IB Interface Bus. A NPS/TPL HG-78K Scanivalve Controller 
with two 48 port Scanivalves allowed the programmed acqui- 
sition of probe and blade surface pressure measurements. 

The software used in the present study for acquisition, 

we reduction and plotting of data were developed from software 

originally created by Duval and Cina. The programs are 
listed and described separately in Ref. 12. 

The uncertainties in the measurements are listed in 


Table l. 


24 


we 


III. EXPERIMENTAL PROGRAM AND RESULTS 


A. PROGRAM OF TESTS 

The test program was in three phases. First, in order 
to verify the new inlet guide vane assembly, tests were con- 
ducted with no bladino in the test section and with the 
upper and lower endwalls set parallel at 35° (design condi- 
tion), 30° and 50° with respect to axial. 

Secondly, tests were made of the effect of wire gauze 
screen materials in reducing non-uniformities in the flow 
entering the test section. Appendix B describes the type 
of screens used and how they were installed. 

The last phase of the test program was a continuation 


of the work initiated by Cina. Table II lists the cascade 


configuration tested. One test was completed successfully 


before aero-mechanical problems were encountered and testing 


was halted until the causes were analyzed. 


B. TEST PROCEDURES 
1. Cascade Adjustments 
In the first and second phases of testing, the same 
procedures were used to realign the cascade for each new 
configuration. The lower and upper end walls were set to 
the desired flow angle and the inlet guide vanes were set 


so that their trailing edges were approximately aligned with 
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the end walls. The flow was started, and the desired inlet 
dynamic pressure was set. All tests were run at an average 
dimensionless inlet velocity (X) of about .13, corresponding 
to an inlet flow dynamic pressure of 18 inches water. Be- 
fore recording data, the water manometer board was checked 
to ensure that the distributions of wall static pressures 

at the inlet plane and outlet plane were acceptably uniform. 
If required, the inlet guide vanes were adjusted to obtain 
uniform static pressure to within +0.5 inches of water. 

In the third phase of testing, initially the pro- 
cedures used by Cina [Ref. 7] were followed, namely: the 
lower end walls were set to the desired inlet air angle and 
the upper end walls were set approximately to the expected 

we exit air angle. The inlet guide vanes were set very approxi- 
mately and the cascade was turned on and set to an inlet dy- 
namic pressure of 18 inches water. The upper end walls and 
the inlet guide vanes were adjusted in turn to obtain wall 
static pressure distributions upstream and downstream which 


were acceptably uniform. Using this procedure however it 


was found on occasion that the inlet air angle sensed by the 

probe at the lower plane at mid-span could be 2 or 3 degrees 
different from the setting of the end walls. 

= The following procedures was subsequently adopted. 

The lower end walls were set to the desired inlet air angle. 

The upper end walls were adjusted to be "wide open", to form 


a diverging passage in which, when the cascade was turned 
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on (to an inlet dynamic pressure of 18 inches water), the 
flow was completely separated. The inlet guide vanes were 
adjusted to obtain the required inlet air angle on the chan- 
nel center line over the center 24 inches in the blade-to- 
blade direction. The upper end walls were then moved 
individually towards the vertical until the lower plane 
static pressure distribution was uniform and the upper plane 
static pressure distribution was acceptably uniform at a 
value close to atmospheric pressure. No readjustment of the 
inlet guide vanes was made. 
2. Measurements 

Probe surveys were carried out in the blade-to-blade 
direction at midspan at the upper and lower planes. In the 
first and second phases of testing, data were taken over 
approximately 24 inches of the test section at intervals of 
0.25 inches. Also, in order to test the repeatability of 
measurements, repetitive samples were taken with the probe 
held fixed at midspan at the lower plane at the center, 10 
inches to the right and 10 inches to the left of center. 

During the third phase of testing, data were taken 


using the procedures established by Cina in Ref. 7. 


C. VERIFICATION OF INLET GUIDE VANE (IGV) ASSEMBLY 
The results of the first phase are presented (as shown 
in Table III) in Figs. 11 to 32. The results are arranged 


into groups. The first group (Figs. 11 to 14) are 
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measurements of tunnel conditions with the end walls set at 
35 degrees. Plots of conditions at the lower plane are fol~ 
lowed by plots of conditions at the upper plane. 

Results for a wall angle of 30 degrees are given next 
(Fig. 15 to Pig. 18), followed by results for a wall angle 
of S50 degrees (Fig. 19 to Fig. 21). Data were taken over 24 
inches at the lower plane, and also at the upper plane at 
30°. At 50°, at the upper plane, only the center 12 inches 
were surveyed. 

The degree of repeatability of conditions in the wind 
tunnel from test to test (with no change in wall setting) is 
demonstrated by the results plotted in Figs. 22, 23 and 24. 

The last group of plots, Figs. 25 to 33, shows the de- 
gree of repeatability in the probe data from scan to scan. 
Data for these plots were obtained by holding the probe sta- 
tionary at midspan in three specific blade-to-blade loca- 
tions in turn and taking 50 repetitive scans of the channels 
normally recorded for survey profile data. The time inter- 


val for each scan was approximately 20 seconds. 


D. TESTING WITH WIRE GAUZE SCREENS 

The selection and installation of the wire gauze screens 
is described in Appendix B. 

The measurements obtained with the various screen con- 
figurations are given (as shown in Table IV) in Figs. 34 to 


45. The results are arranged in four groups. 
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The first group of plots (Fig. 34 to 40) give data ob- 


tained with the 16 mesh .0105 inch diameter wire screen in- 
stalled. Over a blade-to-blade distance from -1.0 to 7.0 
inches (Fig. 35) a peak-to-peak variation in velocity of 
about 1 percent was noted. This is slightly greater than 
the less than 1 percent (0.9 percent) variation noted over 
the same survey region without a screen installed (Fig. 12). 
The results shown plotted in Figs. 41 and 42 were ob- 


tained with two screens installed. One screen (16 mesh, 


-0105 inch diameter) was installed as discussed in Appendix 
B, while the second screen (2 mesh, .0400 inch diameter) was 
attached across the duct at the leading edges of the inlet 
guide vanes. 

The fourth group of plots (Figs. 43 to 46) show the re- 
sults for two different single screens. Probe survey data 
for these screens was taken only at the lower plane. The 
results shown plotted in Figs. 43 and 44 are data obtained 
with a 4 mesh, .041 inch diameter wire screen installed. 

The variation in velocity in the blade-to-blade direction 

was as much as +1.1 percent, veak-to-peak. The results 

shown in Figs. 45 and 46 are for a 5 mesh, .041 inch diameter 
wire screen. The variation in flow velocity was approxi- 


mately +1.5 percent, peak-to-peak. 
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E. PRELIMINARY TESTING OF DCA BLADES 


The results contained in Tables V to IX and Figs. 47 to 


60 are arranged in the following manner. 

The results shown plotted in Figs. 47 to 60 are divided 
into two separate groups. The first group (Figs. 47 to 57) 
contain results which exhibit the quality of the wind tunnel 
flow conditions. The second group (Figs. 58 to 60) shows 
the blade forces (and surface pressures) from survey data. 
In the first group of figures, results are presented first 
to examine the inlet flow uniformity (Figs. 47 and 48); 
second, to examine the outlet flow periodicity (Figs. 49 to 
53); and, finally, to examine outlet flow two dimensionality 
(Figs. 54 to 57). 


All points are shown connected with straight lines. 


30 


Sa i Fe bes = 


we 


IV. DISCUSSION OF EXPERIMENTAL RESULTS 


A. EFFECT OF INLET GUIDE VANE (IGV) MODIFICATION 

The probe survey data shown in Figs. 11 and 12 were 
taken at the lower plane, with the end walls set at 35 de- 
grees. A turning angle of 35 degrees corresponded to the 
"design point" of the inlet guide vanes, when the airflow 
from the plenum chamber was at zero angle of incidence to 
the leading edge of the IGV's. These two figures show that 
the inlet plane total pressure at midspan in the blade-to- 
blade direction had a peak-to-peak periodic variation of 
about +2 percent and therefore about a +l percent peak-to- 
peak variation in the velocity. Corresponding data from 
probe surveys at the upper plane (Figs. 13 and 14) show that 
periodic variations in total pressure were reduced to about 
25% of the value at the lower plane by the mixing of the in- 
let guide vane wakes. 

It can be seen in Figs. 15 to 21 that at "off-design" 
conditions for the IGV's (endwalls at 30° and 50°) there is 
a greater periodic variation in total pressure at the lower 
plane in the blade-to-blade direction than at the design 
point conditions. In Fig. 15 and Fig. 16, it can be seen 
that the periodic variations in total pressure are more pro- 
nounced with the flow from the vlenum at a negative incidence 


angle to the IGV's (endwalls at 30°). Except for the first 
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8 data points in Fig. 15, there is a well defined period of 


about 1 inch of travel. 

Figures 19 and 20 show the probe surveys conducted with 
endwalls at 50 degrees. At this off design condition the 
periodic variation in total pressure was considerably 
greater than the design point and the individual wakes from 
the inlet guide vanes were much less well defined. 

The repeatability of the survey was examined at wall 
angles of both 50° and 30°. Figs. 22 to 24 show that the 
non-uniformities in the flow conditions were repeated to 
(generally) better than 0.5 percent of total pressure. The 
question was then, to what accuracy could the individual data 
points be repeated in successive samples. This was examined 
at several probe positions and the results given in Figs. 25 


to 33 explain the departures in Figs. 22 to 24. 


B. EFFECT OF WIRE GAUZE SCREENS 

All testing with screens was conducted with the end walls 
and inlet guide vanes set to yield a flow angle of 35 de- 
grees. The data obtained with screens installed were there- 
fore compared with the data obtained without screens, shown 
in Figs. 11-14. The effect of the pressure drop across the 
screen on the pressure coefficient plotted in Figs. 34-35 
should be noted. With the first screen installed the drop 
in total pressure from plenum to the probe in the lower plane 


was about 10 inches of water (plenum pressure minus total 
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pressure measured by the probe at the lower plane). With- 
out the screen (at design conditions), the pressure drop 
from the plenum to the probe at the lower plane was approxi- 
mately 2.0 inches of water. Since Q-ef was defined as the 
difference between plenum pressure and lower wall static 


pressure, the value of Q with the first screen installed 


ref 
was about 28 inches water and without the screen installed, 
about 20 inches water. In comparing the peak-to-peak varia- 
tion in Pl seen in Fig. 11 with that obtained with the first 
screen installed in Fig. 34, the difference in the values of 
Qpef must be considered. Examinations showed that the peak- 
to-peak variation in velocity remained at approximately 1% 
when the screen was installed. 

Figures 38 to 40 are plots of data obtained during a 
spanwise traverse of the probe at the lower plane. These 
figures show that the pressure drop through the screen and 
turning vanes was nearly uniform over approximately 8.0 
inches of the 10.0 inch span of the tunnel. 

In an attempt to generate upstream disturbances that 
might trigger early boundary layer transition on the IGV's 
and increase the rate of mixing of the wakes, a second screen 
was attached to the leading edge of the IGV's. The results 
in Figs. 41 and 42 showed that this was not the case and in 
fact the second screen increased the magnitude of the non- 
uniformities at the lower survey plane. Measurements made 


with a single 4 mesh (Figs. 43 and 44) and a single 5 mesh 
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screen (Figs. 45 and 46) of similar blockage showed that 
neither screen influenced the flow in a particularly favora- 
ble manner. The 4 mesh screen caused the peak-to-peak 
variation in velocity to be about +1.1 percent, while the 
5 mesh screen caused the variation to be about 11.5 per- 
cent. This compared unfavorably with the variation ob- 
tained without any screen installed which was less than +1 
percent. 

It was therefore decided to proceed with measurements 


of the test cascade without using screens. 


C. PRELIMINARY TESTING OF DCA BLADES 
1. Inlet Uniformity 
The probe survey at the lower plane in Fig. 48 shows 
that the inlet plane total pressure at midspan varied in the 
blade to blade direction less than 0.5 inches of water, with 
no well-defined spatial period. This was an improvement in 
the inlet conditions found by Cina [Ref. 7: Fig. 16]. That 
the spatial period was not well defined agreed with the 
findings presented earlier in this report. The wall static 
pressure distribution (Fig. 47) showed small variations 
(less than 0.5 ins. water peak-to-peak at the lower plane, 
-4 ins. water peak-to-peak at the upper plane). 
2. Two-Dimensionality 
The data in Figs. 54 to 57 show that, at the down- 


stream plane, an area of (spanwise) nearly uniform conditions 
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existed near the centerline of the cascade. Reference 2 
points out that at higher loadings it is difficult to estab- 
lish a substantial spanwise area of uniform flow in the 
region near the suction side of the blade. This difficulty 
is evident in the data shown in Figs. 54 and 55 which show 
that only about 20% of the spanwise distance is acceptably 
uniform. It is noted that Cina also found reduced areas of 
uniform flow at this incidence angle; however, 30-40% of the 
spanwise distance was found to be acceptably uniform in his 
case. The difference could be the result of the reduced 
spacing of the IGV's and its effect on the side wall boundary 
layers. 

Figure 58 shows results for inlet and outlet flow | 
angles and blade force vectors derived in two ways as shown 
in Appendix B of Ref. 7. These two methods are first the 
applications of momentum conservation to probe survey data 
and second, the integration of surface pressures measured 
over the blade area. Reference 2 points out that for truly 
two dimensional flow the blade forces derived from the two 
methods should be the same. As shown in Fig. 58 the magni- 
tudes and the directions of the two vectors representing the 
blade forces are in reasonable agreement. It is noted how- 
ever that at this particular incidence angle Cina [Ref. 7] 
measured blade forces that were in total agreement in direc- 


tion but disagreed slightly in magnitude. The values of the 
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force magnitudes were about 1.5% lower than those measured 
by Cina. 
3. Periodicity 

As can be seen in Figs. 50 and 51 the total pressure 
and velocity qualitatively repeated fairly well over three 
central blade passages. Acceptably small quantitative dif- 
ferences are noted. There was also a small but measurable 
difference in the surface pressures on adjacent blades, as 
is evidenced in Fig. 49. 

4. Blade Performance 

Figures 59 and 60 are plots of the pressure and 
velocity distributions respectively over the centermost 
blade. These results compare favorably to those obtained 
by Cina for an incidence angle of Pe hee 

Table IX contains the blade performance parameters 
deduced from the probe survey data listed together with the 
data obtained for corresponding test parameters in Ref. 7. 
While differences in two sets of data are evident, the dif- 
ferences are not large. It is noted that the value of the 
loss coefficient was only 10% lower than was measured by 
Cina, but the AVDR was less than 2% different from unity 
rather than the 6.5% measured by Cina. Further measurements 
need to be made, particularly in the light of the following 


discussion, before stronger conclusions can be drawn. 
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5. Aero-Mechanical Problems Encountered 

Fifteen cascade tests were made while evaluating the 
new inlet guide vane assembly and testing the wire screens. 
All runs were made without test blades installed, with a 
plenum total pressure of about 20 inches of water and with 
the upper and lower end walls parallel. No difficulties 
were encountered in establishing the desired flow conditions 
or in using the inlet guide vanes to arrive at a satisfactory 
distribution of wall static pressures. 

The first time the Cascade Wind Tunnel was set up 
with test blades installed to take data at an air inlet angle 
of 39.2°, the tunnel operated normally and the test was com- 
pleted. (The data from this test were subsequently found to 
be highly suspect and are not reported here.) During the 
next test, with an inlet air angle of 42.4, the start-up 
appeared normal and previously established procedures were 
used to arrive at a satisfactorily uniform wall static pres- 
sure distribution. Tunnel operation appeared to be normal 
while taking data, but on shutdown a very noticeable high 
frequency vibration was encountered. Examination of the in- 
let guide vanes revealed that about 40% of the 60 blades 


were damaged. Damage included chips missing from the trail- 


ing edges, blades bent, cracks at the weld where the blade 
is joined to its support and indications that the suction 


side near the leading edge of one blade had been vibrating 


against the pressure side near the leading edge of an ad- 


jacent blade. 


After the inlet guide vanes had been repaired and 
reinstalled extreme care was used at the beginning of the 
next test to adjust the inlet guide vanes, wth two indi- 
viduals monitoring the movement of the adjustment mechanism. 
(IGV adjustment mechanism is described in detail in Appendix 
A.) A lack of stiffness in the mechanism was suspected as 
having been a contributing factor to the failure. 

One successful test was completed at an inlet flow 
angle of 42.4° and these data were discussed above. 

With the next cascade configuration set, at an air 
inlet angle of 45.9°, when the IGV's were adjusted after 

We starting up, high frequency vibrations were again experi- 
| enced. The wind tunnel was shut down and no further testing 
was attempted at plenum total pressures as high as 20 inches 
of water gauge. 

The difficulty encountered with the IGV's is not 
fully understood, however the lack of stiffness present in 
the actuation of the two separate rows of vanes is suspected 
of having allowed the problem to occur. Certainly, the pos- 
sibility of an aerodynamic flutter condition being present 
(due to the misadjustment perhaps) can not be ignored. It 
was noticed after the initial failure thet the lead screw 
which adjusts the IGV's could be turned but the blades 


mounted from only one side would be caused to rotate. This 


could lead to the trailing edge of one blade contacting the 
trailing edge of an adjacent blade and effectively closing 
the blade passages. 

Also, the holes in which the cylindrical shanks of 
the IGV's were held were found not to be uniformly machined. 
As much as 0.1 inches of movement at the tip of some vanes 
was possible while others could barely move. (The most 
seriously damaged blades were found in or adjacent to che | 
larger holes.) 

The tendency for the mechanism to “hang-up" on one 
side would be greater as the vanes became more highly loaded. 
It is noted that the IGV problem was encountered first when 
going to increased incidence angles with the compressor test 
cascade installed. In setting a constant plenum total pres- 
sure of 20 inches of water gauge, the static pressure in- 
crease across the test blades to a constant atmospheric 
pressure at the downstream side implies that a progressively 
increasing dynamic pressure was being generated out of the 
turning vanes. This can be seen in Table IX, where Q) for 


the test at By = 42.4 was 25 inches of water. 


39 


Ce et ae 


V. COMPUTATIONAL PROGRAM 


A. DESCRIPTION OF QSONIC 

The computational code, QSONIC, was developed by the 
staff at the NASA Lewis Research center .? This code is able 
to calculate the blade-to-blade flow conditions in turbo- 
machinery blade rows assuming inviscid flow but including 
streamtube convergence and radius change in the throughflow 
direction. QSONIC is flexible enough to allow the input of 
the appropriate boundary conditions to calculate the flow 
through the test blading in the Subsonic Cascade Wind Tun- 
nel. The program uses a fully conservative solution of the 
full potential equation combined with the finite volume 
method on a body-fitted mesh. QSONIC uses an artificial 
density imposed in the transonic region, if such a region 
exists, to ensure stability and the capture of shock waves. | 

The analysis used by QSONIC is a combination of transonic 
analysis methods to calculate the flow conditions in the vi- 
cinity of a cascade of airfoils. A conservative form of the 
full potential equation is discretized at every point of a 
body fitted periodic mesh and a mass balance is calculated 


through the finite volume surrounding the point. The volume 


lithe help and advice received from Charles Farrell at 
NASA Lewis R.C. in the process of adapting QSONIC to the 
NPS computer is gratefully acknowledged. 
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is corrected three dimensionally for any change in stream- 
tube thickness along a streamtube, if a quasi-3D solution is 
desired. Either elliptic or hyperbolic non-linear partial 
differential equations are used, depending on the local Mach 
number. 
The analysis used in developing QSONIC made the following 
assumptions: 
1) The airflow is inviscid and adiabatic. 
2) The airflow relative to the test blades is steady. 
3) Air is a perfect gas with constant specific heat. 
4) The airflow is isentropic and any discontinuities such 
as shocks are so weak that they may be aporoximated as 
isentropic jumps. 


5) There is no velocity component normal to the 
streamsurface. 


6) The airflow relative to a fixed reference frame (i.e. 
absolute velocity) is completely irrotational. 


Assumption 4 requires that the peak local eiakive Mach num- 
ber on a blade surface be 1.4 or less. The Mach numbers 
measured in test blades in the Subsonic Cascade Wind Tunnel 
would be well within this limit. However, this limitation 
would probably preclude the use of QSONIC for analysis of 
the flow field in the NPS transonic cascade wind tunnel. 
There are some combinations of blading geometry and flow 
conditions which cause unsatisfactory results to be gener- 
ated. For example, because of assumptions 1 and 6, sharp 


leading edges at high incidence angles (more that a few 


degrees) cause large velocity peaks in the blade surface as 
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the flow tries to turn from the stagnation point to the 
suction surface. 

Reference 9 gives a detailed description of QSONIC and 
the solution method used including the governing equations. 
Appendix D describes the operating procedures to use QSONIC 
on the Naval Postgraduate School's IBM 370/3033 computer. 
Appendix D also describes the input and output required as 


applicable to the Subsonic Cascade Wind Tunnel. 


B. APPLICATION TO THE TEST CASCADE 

Appendix D describes in detail the generation of the 
input required for QSONIC when applied to test blading in 
the Subsonic Cascade Wind Tunnel facility. In the present 
work one comparison of code calculations and measured data 
was made before testing was stopped. The comparison was 
for an inlet flow angle (8,) of 42.4°. . 

Tables D.1 and D.3 show the input data generated. Table 
D.6 shows the flow solution output by QSONIC. The flow cal- 
culated on the blade surface, using a 15 by 97 mesh, was 
examined. Figure 61 is a plot of the calculated Mach number 
along the blade surface using two dimensional inputs. Fig- 
ure 62 is a plot of computed Mach number incorporating quasi- 
three dimensional effects. The method of! incorporating 


quasi-three dimensional effects is explained in Appendix D. 
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C. COMPARISON OF CODE CALCULATIONS AND MEASURED DATA 

Table VII lists the data measured in the cascade wind 
tunnel. Co! Cr2 and Xvel are defined in Appendix C. The 
surface Mach number distribution measured on the center blade 
is shown plotted in Fig. 63. 

For comparison, the computed two dimensional, comouted 
quasi-three dimensional, and the Mach number measured in the 
cascade wind tunnel are plotted together in Fig. 64. 

Excellent agreement between all three cases is seen. 

As would be expected, the greatest difference between mea- 
sured and calculated data is near the leading edge in the 


suction side and at the trailing edge of the blade. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


Based on the results of the first part of the present 


study, to evaluate the effects of the altered inlet cuide 


vane spacing on flow uniformity and periodicity, the follow- 


ing conclusions were drawn: 


1. 


With the inlet guide vanes operating at design, the 
peak-to-peak variation in velocity was about 11%, and 
there was a well defined spatial period of about 1 
inch. 


Operating the cascade wind tunnel in a configuration 
that requires the inlet guide vanes to be set to other 
than zero incidence resulted in peak-to-peak varia- 
tions in velocity greater than +t1% and a spatial period 
that was less well defined. 


The second part of the study, to evaluate the use of 


wire gauze screens to further reduce the non-uniformities in 


the flow field, led to the following conclusions: 


1. 


A 16 mesh screen with a blockage factor of .69 had a 
slightly aggravating effect in the variation in velo- 
city. The peak-to-peak variation in velocity with the 
screen installed was slightly greater than with no 
screen installed. This occurred at the expense of a 
pressure drop of about 10" of water across the screen. 


The use of screens with similar blockage but with 
larger mesh and larger diameter wires resulted in 
larger peak-to-peak variations in the velocity at the 
lower plane. 


The overall objective of the present study was to mea- 


sure the performance of the DCA test blading. Because of 


aero-mechanical problems encountered with the inlet turning 


vanes the performance of the blades was obtained at one 
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incidence angle only. The following was concluded from the 


limited test program: 


1. 


we 


As a result of the reduced inlet guide vane spacing 
the variations in velocity and total pressure at the 
inlet plane were much less than those reported by Cina. 


Good periodicity was found from one blade passage to 
another. 


An acceptable region of spanwise uniformity (20-40% of 
blade span) was found at the downstream plane at the 
one test condition reported. However, this was less 
than was previously reported for the same incidence 
angle. 


The blade forces derived from the integration of sur- 
face pressure measurements and probe survey data were 
in close agreement in both magnitude and direction. 


The Mach number measured by surface pressure taps over 
the surface of the blade and Mach number calculated 
using the program QSONIC were in excellent agreement 
qualitatively and reasonable agreement quantitatively. 


The specific reasons for the aero-mechanical problems 
experienced with the inlet turning vanes have not been 
identified completely. 


Based on these conclusions and other observations, the fol- 


lowing recommendations are made: 


p 


Use of the present inlet guide vane assembly and ad- 
justment mechanism for testing at inlet dynamic pres- 
sures higher than about 15 inches of water is unsafe. 
There are three possible solutions to this problen. 


a) Operate only within the dynamic pressure range of 
10-15 inches of water. 


b) Modify the new inlet guide vane assembly so that 
the vanes are supported at both ends on their axes 
of rotation. [Supporting the IGV's from both ends 
would prevent flapping vibrations of the (present- 
ly) cantilevered vanes. Such vibrations, when the 
vanes are supported alternately from opposite ends 
and the gaps are small compared to the chord, 
might lead to a potentially destructive flutter 
mode at particular flow velocities.] 
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c) Replace the entire inlet turning vane section with 
one of entirely new design. 


The procedure should be adopted immediately of adjust- 
ing the vanes and walls of the cascade at lower values 
of the dynamic pressure before increasing the blower 
speed to the desired operating condition. 


More time needs to be spent to examine the flow field 
produced between the guide vanes and the test blades, 
and to establish the effects of the movement of the 
tail boards. The uniqueness of the flow field when the 
wall static pressure uniformity is used as a criterion 
of good inlet flow, needs to be examined by conducting 
repetitive tests at nominally similar test parameters. 
Only when the adjustment of the flow in the facility 
and the quality of the flow itself is fully understood 
should the measured blade performance data be accepted 
as final. 


The upper electrical yaw adjustment mechanism should 
be replaced with a manual system to greatly decrease 
the time required to achieve probe pressure (angle) 
balancing. 


Develop the computer code necessary to take advantage 
of the plotting data created by QSONIC. 


Modify the data acquisition and reduction software for 


the HP 9845 so that real time plots of blade perfor- 
Mance parameters can be displayed. 
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TABLE I. MEASUREMENT UNCERTAINTY 


Item Description Method Uncertainty 
x Blade-to-Blade dimension Position +.01 in. 
x = Oin. West end Potentiometer 


x = 60 in. East end 


Zz Spanwise dimension Position +.01 in. 
z= Oin. North wall Potentiometer 
z= 10 in. South wall on probe mount 
By Inlet flow yaw angle Angle Potenti- +.2 deg. 
ometer on probe 
mount (hand 
adjustment) 
8. Outlet flow yaw angle Angle Potenti~ +.5 deg. 
ometer on probe 
mount (motor 
driven adjustment) 
P len Plenum total pressure Static tap in +.01 in. H,0 
P plenum chamber gauge 
vV=0 
P Static pressure at the Calibrated pneu- +.1 in. HO 
test plane matic probe gauge 
Pie Static pressure at Static tap on +.Q1 in. H,0 
x = 0 in., y = -16.25 in., North wall gauge 
z= 9 in. 
P ATM Atmospheric pressure Absolute Strain +.3 in. H,0 
Gauge Transducer 
P Pressure Scanivalve +.01 in. H,0 
Transducer gauge 
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TABLE II. CASCADE CONFIGURATION FOR DCA BLADE TESTS 


Constant Parameters 


Number of Blades 20 
Spacing (Pitch) 3 inches 
Chord 5.01 inches 
Solidity 1.67 
Thickness 7.0 percent of chord 
ve Camber Angle 45.72 degrees 
- Stagger Angle 14.72 degrees 


Variable Parameters 


42.4 degrees 


i 5.3 degrees 
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By Survey Plane Survey Direction Fig. Nos. 


35 


30 


50 


50 


30 


30 


TABLE III. 


Lower 


Upper 


Lower 


Upper 


Lower 


Upper 


Lower 


Lower 


Lower 


SUMMARY OF MEASUREMENTS WITHOUT SCREENS 


(24 


(24 


(24 


(24 


(24 


(12 


(24 


(24 


B-B 
inches) 


B-B 
inches) 
B-B 
inches) 
B=-B 
inches) 
B=-B 
inches) 
B=-B 


inches) 


B-8 
inches) 


B-B 
inches) 


Fixed Probe 


(10" 


L. of ¢) 


(on €) 


(10" 


R. of €) 
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ll & 12 


13 & 14 


15 & 16 


17 & 18 


19 & 20 


21 


22 & 23 


24 


25-27 


28-30 
31-33 


Purpose 


Flow Field 
Determination 


Survey 
Repeatability 


Point 
Repeatability 
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TABLE IV. SUMMARY OF MEASUREMENTS WITH SCREENS 


=o ro) 
(By = 35°) 
Screen Survey Plane 
16 mesh Lower 
-0105 wire 
Upper 
Upper 
16 mesh Lower 
-0105 wire 
+ 2 mesh ‘ 
ahead of IGV's 
4 mesh Lower 
-041 wire 
5 mesh Lower 
.041 wire 
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Survey 


Direction 


Spanwise 


B-B 


B-B 


Fig. Nos. 


34 & 35 
36 & 37 
38-40 


41 & 42 


43 & 44 


45 & 46 


TABLE V. PROBE DATA, UPPER PLANE AT MIDSPAN (i = 5.3°) 


wnTA FROM File WEREDe: Tis 
BLADE TO BLADE TRAVERSE MLODPAH 


WRPFER PLANE 


Fotme Locsin: Beta - Qlbar Fs-lbar Fr mtbar re Me ae 
ROKER ECRRSEK THREE CEBSEREERERREE RRR SEER RRRAREERER EKER EEO SEH RE KOH 
-7.31 “1.59 » 5686 
-6.34 -.34 »3532 
-6.3? -.33 04276 
-3.3? -2.56 5696 
“5.67 ~2.56 Sez 
-5.43 -2.56 ~57°36 
-S.23 -1.79 -5739 
-S.a4 -1.79 »5747 
-4.35¢0 -t.78 - 35697 
-4.66 -1.7%3 -5785 
-4.4t -1.72 5687 
“4.15 -1,36 - 5658 
“3.96 -1.35 - 5629 
-3.73 -1.38 5462 
-3.53 -1.36 -4863 
3.39 -1.35 - .4349 
19 -1.88 24239 
-3.16 0 4836 

-2.41 .S46?7 

-2.193 3717 

-1.92 | .5751 

«5773 

S814 

5312 

~5765 

goede 

25748 

3679 

» 5378 

04.704 

24148 

24295 

*.5611 

S362 

oor76 

» 5795 

- 5880 

3915 

-5902 

» 3983 

- 3816 

24499 

4421 

~S72e2 2688 

- 28763 0259" 
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ws 
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mG Ge fh 
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- 
Sw ery oh & oF fy 
Wa od fe GA 
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PMR OAS, 
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AAW AOA Mh | 


Ot ay 


HBOwWOGRKe Aaa & Meo Ga G os he 
PQ HDMK Oo wIn Desk HAS w te @ 


3 


29 
21 
22 
23 
24 
25 
26 
2? 
23 
29 


AME EELARNEEEWOAAGU EEE ELE 


70, 


TM oA, 5 
SD IAMAMMACHAH EES 


cone hk enmmnk. SON 
ee fees Tm ON & = OG Or ue 


a, 


TABLE VI. 


BLADE Ta BLAULE TRAVERSE MIDSPAN 


LOWER PLANE 


PROBE DATA, LOWER PLANE AT MIDSPAN (i 
SATA FROM FILE LRREDSST14 


5.3°) 


Point Loctind Beta He Gitar Pe“Olbar Fregloar . bar 

RSL RREREK KEKE REE ARE EH SCKREREREREELERERREKSREERERERERE CRRA RS RRO 
1 -4.48 ~42.45 239533 -.,1372 8932 345? 
< -3.56 “42.42 ~ 9603 -. 1444 8967 3474 
3 -3.588 -42.42 -3567 -.1427 4936 2454 
4 -2.58 ~-42.44 «9614 -.1519 1631 . 3488 
5 -2.08 $42.43 09741 -.14692 agri .3524 
5 -1,50 -$2.43 .39767 -.1493 6375 23532 
7 -1.88 -42.43 ~9722 -.1481 8883 .3522 
a -.58 -42.43 .9743 -.1364 3921 » 3546 
3 0.98 -42.43 29799 -. 1567 6935 .8549 
ig -3@ -42,.42 - 3849 -.1592 ~oss4 .3585 
11 1.09 -42.43 ~ 9814. -. 1684 6 6858 .3593 
le 1.56 -42.43 9671 -.1628 8779 »8635 
13 2.88 -42.43 97387 -.1613 . 8898 e582 
14 2.56 -42,43 29742 -.1632 2971 .35355 
15 3.028 -32.41 973i -.16358 . 1059 .8526 
16 3.58 -42.42 ‘988? -.1663 ~ 83924 »33589 
1? 4.00 -42.44 ©9738 -.1631 . 8968 »356? 

“4 
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TABLE VII. 


pe 
eo + 
& 


0 wo 


Coe ER > fe 
Ofte Uh 
Ce Oe co 
Grew we 


orn? «) 


vee Cpl 


PRESSURE SIDE CENTER BLADE 


6os4 -6655 
Beals ~5417 
BBES 3240 
Bitz 23287 
9215 24871 
- 9383 24813 
. 3452 24768 
-6576 » 4737 
- 8663 ~ 4764 
.8715 24892 
8736 «4871 
8727 04771 
4073 24956 
3601 - 4839 
0487 24949 
0411 4672 
6327 4537 
382328 24245 
BLs3 ~3815 
8aae 22778 


SULTION SITE CENTER BLADE 


.022? =1. 3641 
.8210 -. 7756 
.9389 -. 4986 
. 6563 -.3541 
.8719 -.3491 
.3978 -.2923 
.1176 =. 2860 
1309 -. 2038 
11399 -.19540 
.1432 ~ 20854 
J1412 -.0584 
£1339 . 8102 

.1269 .8631 
.18Z1 1456 
8695 219i 
.9755 2326 
.0593 2636 
9407 . 2842 
£8286 .2948 
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TABLE VIII. ADJACENT BLADES DATA (i = 5.3°) 


wee Wee Tpl Cpe Mach atte | 


RTRE OK FRRES KERESR ERE PS LRRERREARA HRS REEREEESHRERRREEEE ROR RARE 


PRESSURE SIDE LEFT BLADE 
£1218 9383 64014 -.0251 2462 1068 
4192 18716 4619 0568 2275 1012 
3283 8414 4491 8635 . 2302 1024 
SUCTION SIDE LEFT BLADE 
i219) « 0710 -.3438 “1.3791 3660 1615 
4192 1468 -.1505 -1.8273 3371 1491 
3283 8895 1975 -.3943 2794 1246 
PRESSURE SIDE RIGHT BLADE 
1218 2303 4648 8907 2278 1618 
4192 8716 4658 8939 2266 . 1068 
22383 0411 4597 8329 2279 1014 
SUCTION SIDE RIGHT BLADE 
61218 8718 -.3z85) -1.3513 . 3638 1606 
14132 . 1400 -.1519 -1,8299 3374 11492 
8233 8395 -19L1 -. 4989 . 2885 1245 
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TABLE IX. BLADE PERFORMANCE 


Present Results 


By 42.42 
i 5.3 

6 10.44 
D 0.455 
wo 0.037 
w cos” B. 

——_;- 0.020 
20 cos By 

Ww cos B45 2 

—j4 (x 10°) | 1.09 
AVDR 1.015 
CPsraric net? 
CoM -1.385 
Com -0.669 
Cop -0.572 
Q, (in. H,0) 25 
xX ~l4 
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From Ref. 7 


42.42 
5.3 
0.4 
9.0 
0.46 
0.041 


0.023 


1.242 


1.065 
0.351 
-1.380 
-0.566 
-1.476 
-0.645 
22 
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Fig. 7. Photograph of Instrumented Blade 
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Fig. 8. Instrumented Blade 
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Fig. ll. Probe survey Data at Midspan of Lower Plane 
End Walls at 35°, Points 1 to 50 
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Fig. 13. Probe Survey Data at Midspan of Upper Plane 
End Walls at 35°, Points 1 to 50 
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Probe Survey Data at Midspan of Upper Plane 
End Walls at 35°, Points 51 to 100 
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Fig. 15. Probe Survey Data at Midspan of Lower Plane 
End Walls at 30°, Points 1 to 50 
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Fig. 16. Probe Survey Data at Midspan of Lower Plane 
End Walls at 30°, Points 51 to 100 
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Fig. 17. Probe Survey Data at Midspan of Upper Plane 
End Walls at 30°, Points 1 to 50 
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Fig. 19. Probe Survey Data at Midspan of Lower Plane 
End Walls.at 50°, Points 1 to 50 
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Fig. 21. Probe Survey Data at Midspan of Upper Plane 
End Walls at 50° 
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Fig. 22. Probe Survey Data at Midspan at Lower Plane 
End Walls at 50°, Two Runs, Points 1 to 50 
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Fig. 23. Probe Survey Data at Midspan at Lower Plane 
End Walls at 50°, Two Runs, Points 51 to 100 
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Fig. 25. Repetitive Samples with Fixed Probe Position 
(10" Left of CTR Midspan, End Walls 30°, 


Lower Plane (POL ENUM - Pama) /Qre¢) 
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Fig. 26. Repetitive Samples with Fixed Probe Position 
(10” Left of CTR Midspan, End Walls at 30°, 
Low 2 - 
ower Plane (PY Pama) /Q+e¢) 
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Fig. 27. Repetitive Samples with Fixed Probe Position 
(10" Left of CTR Midspan, End Walls at 30°, 
Lower Plane (PoLENUM - Pi)/Ore¢ 
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Fig. 28. Repetitive Samples with Fixed Probe Position 
(on Centerline at Midspan, End Walls at 30°, 
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Fig. 22. Revetitive Samples with Fixed Probe Position 
(on Centerline at Midspan, End Walls at 30°, 
3 Lower Plane (PY - Pamp) /Qre¢) 
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Fig. 30. Repetitive Samples with Fixed Probe Position 
(on Centerline at Midspan, End Walls at 30°, 
Lower Plane (POL ENUM - Pi)/Qe¢) 


84 


Piplenumi-PC amb) -“QreFf 
14 INCHES PT OF CTR 94 DES LOWER 


er 


i 
’ 
‘ 
‘ 


ee ee ee ee 
We t 
{ 


F 


J | FO Wael a Snes Ce 
19; ae ee Ce OU CS US Rn GOES SESW GOON Sa Sn Om SOOT Mn Se GONDD AOY GORS CY COUSCOUS CT IS DON SD LE CO KOON) UO DN JW Sant 
‘ ws iv) ” a) Yor wow © MBM wr eo Dt & 
wo «= Q@Qow+ = a ry wot nN rs 


3a 


Fig. 31. Repetitive Samples with Fixed Probe Position 
(10" Right of CTR Midspan, End Walis at 30°, 


Lower Plane (PoLENUM - Pamp)/ res) 
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Fig. 33. Repetitive Samples with Fixed Probe Position 
(10" Right of CTR Midspan, End Walls at 30°, 
Lower Plane (Poy ony 7 Pt) /Oree) 
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Fig. 35. Probe Survey Data at Midspan, Lower Plane 
16 Mesh Screen, Walls at 35°, Points 51 to 100 
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Fig. 36. Probe Survey Data at Midspan, Upper Plane 
16 Mesh Screen, Walls at 35°, Points 1 to 50 
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Fig. 38. Probe Survey Data Span Traverse, Lower Plane 
16 Mesh Screen, Walls at 35°, 10" Left of CTR 
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Fig. 40. Probe Survey Data Svan Traverse, Lower Plane 
16 Mesh Screen, Walls at 35°, Center of Test 
Section (POLENUM - Pi)/Q-e¢ 
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Fig. 41. Probe Survey Data at Midspan, Lower Plane 
16 Mesh and 2 Mesh, Walls at 35°, Points 1 to 50 
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Fig. 42. Probe Survey Data at Midspan, Lower Plane 
16 Mesh and 2 Mesh, Walls at 35°, Points 51 to 100 
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Fig. 43. Probe Survey Data at Midspan, Lower Plane 
4 Mesh Screen, Walls at 35°, Points 1 to 50 
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Fig. 44. Probe Survey Data at Midspan, Lower Plane 
4 Mesh Screen, Walls at 35°, Points 51 to 100 
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Fig. 45. Probe Survey Data at Midspan, Lower Plane 
3 Mesh Screen, Walls at 35°, Points 1 to 50 
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Fig. 46. Probe Survey Data at Midspan, Lower Plane 
5 Mesh Screen, Walls at 35°, Points 51 to 100 
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Fig. 47. Wall Static Pressure Distribution 
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Fig. 48. Probe Survey Data at Upstream Midspan 
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Fig. 52. Probe Survey Data at Midspan 
(i 5.35 (P - Py /21° Upper Plane) 


106 


GUTLET ANGLE 
UPPER PLANE MIDSPAN (1-5.35 


1.5 


| 
| 
| 
4 
| 
| 


rd 


ab 


if 


[INCHES 


53. Probe Survey Data at Midspan 
(i = 5.3, Outlet Angle, Upper Plane) 
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Fig. 57. Spanwise Probe Data Surveyed 1 in. from 
“7 Pressure Side of Centermost Blade 
a (i = 5.3, X/X, Upper Plane) 
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Fig. 58. Resultant Blade Force Vectors by Momentum 
Balance ) and from Surface Pressure 
Integration (————) Ls 5.3 
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Fig. 59. Measured Blade Surface Pressure Distribution 
(i = 5.3, * = Pressure Side, 
+ = Suction Side) 
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3D Code Blade Surface Mach Number Distribution 
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APPENDIX A 


MODIFICATION TO THE INLET GUIDE VANE SECTION 
OF THE SUBSONIC CASCADE WIND TUNNEL 


As discusse@ in Section I, Cina discovered during his 
test program that while the inlet flow to the test section 
was uniform in direction and uniform in wall static pres- 
sure, it contained a variation in velocity and stagnation 
pressure resulting from the wakes of the IGV's. Because the 
inlet guide vanes were spaced at intervals of two inches and 
the test cascade blades spaced at three inches, departures 
from strictly periodic conditions were detected from one 
test blade passage to another. 

To alleviate this problem, the inlet guide vane arrange- 
ment was modified so that the guide vanes were placed at l 
inch intervals. In ’order to preserve the option of reverting 
to a two inch IGV arrangement and because it was not possi- 
ble to machine the south wall to hold additional blades, a 
separate structure was placed between the bell mouth contrac- 
tion and the walls of the cascade. By mounting the IGV's 
in a separate unit which remained fixed once installed, 
hardware adjustments between tests associated with a change 
of end wall angle were greatly simplified. 

The new inlet guide vane assembly was constructed using 


two lengths of 10 inch steel channel as shown in Fig. A.1l. 
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One set of guide vanes was mounted on the south side of the 
unit at 2 inch spacings. A second set of vanes was mounted 
at 2 inch spacings on the north side of the unit. When the 
unit was assembled the guide vanes mounted from alternate 
sides meshed, resulting in an inlet guide vane spacing of 

1 inch. The unit was provided with a single hand crank at 
the east end so that the vanes would be adjusted in unison. 
Once installed the complete structure could be left in place 
when the cascade north wall was removed to adjust air inlet 
angle. The one inch vane spacing ensured that periodicity 
at the test section would result for any test blade spacing 
which was a multiple of 1 inch. Equally important, the 
wakes remaining at the inlet to the test cascade would be 
greatly reduced as a result of closer spacing. 

Figure A.1 shows the details of the inlet guide vane 
unit, while Fig. A.2 shows the assembly in relation to the 
bellmouth contraction and the side walls. Figure A.3 shows 
the mechanism to adjust the inlet guide vanes. Figure A.4 
shows a view of the IGV assembly from the north side. A 
view of the Cascade Wind Tunnel partially assembled (north 


side wall off) is shown in Fig. A.5. 
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Fig. A.1. Inlet Guide Vane Assembly 
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APPENDIX B 


SELECTION AND INSTALLATION OF SCREEN MATERIAL 


Pankhurst and Holder [Ref. 17] show that the turbulence 
of an airstream can be increased by placing a coarse mesh 
across the flow upstream of the test section. One of the 
most effective methods which is used for the reduction of 
turbulence and non-uniformities also consists of placing a 
mesh screen across the tunnel. Screens used for this pur- 
pose are of a finer mesh and are placed at a greater dis- 
tance from the test section, and normally in the low-speed 
region upstream of the contraction in a conventional sub- 
sonic wind tunnel. By using such a screen the large scale 
eddies are removed at the expense of the introduction of a 
greater number of smaller eddies which decay rapidly. 

McEligot [Ref. 16] investigated the possibilities of 
reducing non-uniformities in the test section of the sub- 
sonic cascade wind tunnel. His investigation and recom- 
mendations were completed while the inlet guide vanes were 
still at 2 inch spacings. McEligot concluded that some 
modification was necessary to achieve one percent uniformity 
for the mean velocity at the test cascade inlet plane and 
suggested several options. One of the options suggested was 
placing the turning vanes (inlet guide vanes) at a closer 


pitch. As explained in Appendix A, the pitch of the inlet 
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guide vanes was reduced from 2 inches to 1 inch. This new 
inlet guide vane arrangement did result in a one percent 
uniformity for the mean velocity at the test cascade inlet 
plane. 

The other approach suggested was the use of wire gauze 
screens. McEligot showed that the velocity distribution ap- 
peared to be largely dependent on a pressure drop coefficient 


K, defined by the equation 


P, ~ P2 
aa crepe, [a 


OV 
where Py and P» are the pressures upstream and downstream of 
the screen respectively. This pressure drop coefficient de- 
we pends mainly on the blockage coefficient 8 defined by the 


equation 
B = (1 - az)? 


where d is the diameter of wire used in the screen and & is 
the distance between the wires. This blockage coefficient 
is commonly referred to as “percent open area" in catalogues 
of industrial wire cloth and woven wire screens. 


For the velocities and flow angles used in the cascade 


wind tunnel, McEligot recommended using a wire gauze screen 
with a resistance coefficient, K, of 2.2 and a blockage 
coefficient, 8, of 0.47. However, since the new inlet guide 


vane arrangement resulted in a one percent uniformity for 


4th 


ve 


the mean velocity and the pressure drop across a screen with 
a blockage coefficient of 0.47 was expected to be higher 
than could be tolerated for the desired test conditions, 
screens with a slightly higher blockage coefficient (higher 
percent of open area) were selected to be tested. 


The screens tested were of the following configurations: 


MESH WIRE DIAMETER (inches) BLOCKAGE COEFFICIENT 
4 -0410 - 6989 
5 -0410 - 6320 
16 0105 -6922 


Until the effectiveness of wire gauze screen in reducing 
non-uniformities in this facility was proven, a temporary 
installation of the screen material was considered adequate 
for testing purposes. The test screen was installed in the 
cascade wind tunnel by placing it between the inlet guide 
vane assembly and the north and south end walls. This ar- 
rangement placed the screen 7.25 inches downstream of the 
inlet guide vanes and19.3 inches upstream of the lower test 
Plane. Figure B.1 shows the installation of the wire gauze 


screen. 
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APPENDIX C 


CASCADE PERFORMANCE PARAMETERS 
(by F. S. Cina; reproduced with minor changes from Ref. 7) 


The performance of a cascade is specified in terms of 


the deviation angle (6) and the loss coefficient (w) for 


given inlet conditions. In Ref. 1 the loss coefficient is 


shown to correlate in terms of the Diffusion Factor (D). In 


the present work, the performance parameters were calculated 


using the following expressions: 


1. Loss Coefficient (w) 


where the mass averaged pressure coefficients 


defined in Appendix C of Ref. 7. 


C of Ref. 7 that the effect of time dependent 


ditions are removed and the effect of AVDR is 


the use of Eq. (1). 


2. Diffusion Factor (D) 


pie 
W ~ 2oW, 
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It is shown 


(1) 

es 
in Eq. (1) are 
in Appendix 
supply con- 


included in 


(2) 


3. Pressure Rise 


Cc = eset (3) 
Pstatic Qy 


4. Blade Surface Pressure Coefficients 


Pp - BD 
Qy 
Pp - Dp 
C2 7 = (5) 
Qo 


5. Dimensionless Velocity 


Si Se 
YE V3CT, 


where V is the local velocity, Vv. = v2C oT. is the "limiting" 


velocity and tT, is the stagnation temperature. 
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APPENDIX D 


INSTRUCTIONS FOR PREPARING INPUT AND OPERATING 
QSONIC USING A RECTILINEAR CASCADE CONFIGURATION 


D.1 BACKGROUND INFORMATION 

QSONIC has the capability to calculate an axial, mixed 
or radial flow field and the test cascade can be rotating or 
stationary. The geometry of the streamsurface can be a 2D 
planar cascade or axisymmetric with varying channel thick- 
ness and radial position. The capabilities of QSONIC, be- 
yond those of previous cascade analysis methods (such as 
described in Ref. 8) include the ability to calculate 
through weak shocks with a peak relative Mach number less 
than 1.4, and completely around both leading and trailing 
edge regions of a blade profile. The blade shapes in the 
leading and trailing edge regions are not restricted to 
circular arcs. Detailed instructions for preparing input 
for a configuration other than an axial flow, stationary and 
rectilinear cascade may be found in Ref. 9. What follows 
are instructions for preparing the input applicable to the 
Rectilinear Cascade facility and running QSONIC on the Naval 
Postgraduate School's,computer and associated operating sys- 
tem. It is assumed that the reader has a working knowledge 
of the NPS computer operating system and is familiar with 


Refs. 13 and 14. 
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QSONIC operates in two parts. The first part gener- 
ates a body (blade) centered mesh (geometry generation). 

The second part actually solves for the flow conditions at 
points in the mesh (flow solution). The data necessary to 
generate a mesh consists of a two-dimensional description of 
the blade shape. This is in the form of pairs of (X,Y) 
points on the surface, together with parameters that de- 
scribe the cascade layout, such as chord and stagger angle. 
Additionally, parameters describing the density of mesh 
lines complete the input for the geometry generation. 

For flow field calculations, the upstream flow condi- 
tions, convergence criteria and a schedule of meshes to be 
used should be input. If quasi-three dimensional effects 
are to be considered, a data file containing a description 
of the streamchannel's radial thickness and position as a 
function of distance along the stream surface is needed. 

For the case described herein, this was input by assuming a 
linear reduction in streamchannel thickness using a factor 
of 1/AVDR. This gave excellent results. (The output of 
another NASA code, Meridl [Ref. 15], can be used to input 
data to QSONIC for compressor flow field calculations. This 
program has recently become operational on the NPS computer.) 

The output of QSONIC consists of listings which con- 
tain an echo print of the input data, generated mesh coordi- 
nates on the blade surface, progress reports on the flow 


convergence and a list of the final velocities, pressures 
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and densities on the blade surface for each grid that was 


included in the schedule of solution meshes. 


D.2 INPUT DESCRIPTION 

The input for QSONIC falls into the following 
categories: 

Logical and case control parameters (NAMELIST PARAMS) 
Bulk data input: 

For geometry generation runs (NAMELIST INSTUFF) 

For flow solution runs: 

Mesh point storage files 

Streamchannel data file (for quasi-3D) 

The following is a description of the logical and case 
control parameters for the Rectilinear Cascade. Except for 
TITLE, the format for all these variables is in namelist 
form. The namelist is PARAMS. This information is taken 
from Ref. 9 and adapted to the Cascade Wind Tunnel. Name- 
list variables can be entered in any order. If a default 
value is listed, it is not necessary to enter that particu- 
lar variable. If the default value is listed as none, then 
a value for that parameter must be input. 

The following parameters apply to both the mesh gen- 
eration and the flow solution runs, but the values need not 


be the same. 
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Variable 


Name 


TITLE 


NOFLOW 


NOZES 


Type 


Alphameric 


Logical 


Integer 


Integer 


Default 
Value 


None 


False 


Description 


This is a one-line name for 
the case being run. TITLE 
must appear on the first line 
of the data file that will be 
referred to as NAMELIST DATA. 


NOFLOW = .TRUE. if a run is 
to stop after generating a 
mesh, such as the mesh genera 
tion run. 


NOFLOW = .FALSE. for the flow 
solution run. 


MS is an array (max. dimension 
10) of values (mx. value = 
25) for the number of grid 
lines in the mesh that will 
enclose the blade. 25 is a 
Satisfactory number for blades 
of solidity near unity. As 
solidity increases the maxi- 
mum value in MS should de- 
crease. 


NOZES is an array (max. 
dimension 10) of values (max. 
value = 49) for the number of 
grid lines in the mesh ra- 
diating from the blade on one 
surface. If a value of NOZES 
is greater than zero, a mesh 
with that many lines will be 
developed and stored in the 
file MESHGEN DATA. If the 
value is negative, QSONIC 
assumes that this MESH al~ 
ready exists in the file 
MESHGEN DATA and will be read 
in. For Geometry Generation 
runs NOZES > 0 and for flow 
solution runs NOZES < 0. For 
electrostatic analog grid 
generator, NOZES must be odd. 


BETAL Real 


e 
t we BETA2 Real 
® 

e 

GAMMA Real 
@ 

TOLS Real 
e 


Bak eet ee ee en gt Ce ae Cee TOT ee a Fr ae 


None 


None 


1.4 


None 


BETA1 is the flow angle at the 
upstream boundary, in degrees. 
BETA] is measured from the aero- 
dynamic chordline to the direc- 
tion of flow: clockwise is 
negative. This can be obtained 
from the output of the cascade 
tunnel data reduction program 
'CX4431'. It is listed as "in- 
let air angle." Because of a 
difference in definition, it 

is necessary to use the usual 
inlet air angle minus the stag- 
ger angle for BETA1. 


Example: The output of CxX4431 
lists an inlet air angle of 
42.429°. The cascade is con- 
figured with a stagger angle of 
14.27°. 


42.429° 
- 14.27 


28.16 = BETAI1 


BETA2 is the flow angle at the 
downstream boundary in degrees. 
It is measured from the aerc.: 
dynamic chordline to direction 
of flow. Clockwise is negative. 
This can be obtained from the 
output of the cascade tunnel 
data reduction program 'CX4431'. 
It is listed as outlet air 
angle. Because of a difference 
in definition, it is necessary 
to use the outlet air angle 
minus the stagger angle for 
BETA2. 


This is the ratio of specific 
heats. For the Subsonic Cas- 
cade Wind Tunnel the default 
value works well. 


This is an array of dimension 
10 of tolerances corresponding 
to MS and NOZES. Each grid 
solution will proceed until its 
TOLS value is satisfied. 
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MESH1 


MESHN 


LAMDAO 


CHORD 


Integer 


Integer 


Real 


Real 


Real 


None 


yone 


None 


None 


This is an index in the 
arrays MS and NOZES of the 
first mesh to be generated 
and/or used for the flow so- 
lution. For geometry genera- 
tion runs, MESH1 selects one 
of the grids to be stored, 
provided NOZES(MESH1) > 0. 


This is an index in the 
arrays MS and NOZES of last 
grid to be calculated for 
geometry generation runs or 
used for the flow solution. 
For geometry generation runs, 
MESHN = MESH1. Subsequent 
flow solution runs then solve 
the case for all grids listed 
in MS and NOZES between index 
MESH1 and MESHN. For fJ.ow 
solution MESH1 normally 
points to the coarsest mesh. 


To insure equal spacing of 
grid lines over the entire 
mesh, 


where MESH1 < I < MESHN and 
K= 2, 3, 4, . .. (K = 2 if 
grid lines are doubled be-~- 
tween successive grids). 


Stagger angle of the blade 

row in degrees measured from 
the throughflow direction to 
the blade chord line. (Clock- 
wise is negative.) 


True chord of the blade, in 
the same units as the blade 
coordinates. 


Blade spacing in the same 


units as the blade 
coordinates. 
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The following parameters are required only for the 
geometry generation run. QSONIC gives the user a choice of 
two grid (mesh) generators. The Electrostatic Analog grid 
generator is applicable to any blade shape and any value of 
stagger or turning. The interpolation scheme grid genera- 
tor works for most blades except those with sharp leading 
edges where the edge radius is less than .5% of chord. The 
interpolation scheme allows the user to concentrate grid 
lines in areas of high interest around the blade. With the 
Electrostatic Analog grid generator no concentration of grid 
lines is available. For the cascade configuration used in 
this study, the interpolation scheme provided gross errors 
in the flow solutions, so the Electrostatic Analog grid 


generator was used with good results. 


Name Type Default Description 
NED Integer None Total number of body definition 


coordinates that are input. These 
are the X,Y points that describe 
the profile of the blade, with the 
first point repeated as the last 
point. 


KN Integer None KN is used to indicate which grid 
generator is to be used. KN = 0 
will call the electrostatic analog 
grid generator. For interpolation 
scheme, KN = number of body points 
on upper surface from the minimum 
to maximum X points, inclusive. 


NED and KN are required for eit':': grid generator. If 


the electrostatic analog grid generator is used, no other 
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parameters are used. 


If this grid generator is used the 


value from NOZES that is used must be odd. 


The following parameters are used only if the Interpo- 


lation Scheme is used to generate the mesh. 


Name 


RLE 


THETL 


THETT 


CAMPER 


STABAC 


Type Default 
Real None 
Real None 
Real 0 
Real 0 
Integer 6 
Real 0.999 


Description 


RLE is the leading edge radius 
of the blade, with units the 
same as the coordinates de- 
fining the blade profile. 


RTE is the trailing edge radius 
of blade. 


THETL is the camber angle of 
the leading edge in degrees. 

It is measured from the aero- 
dynamic chord to the line tan- 
gent to the mean camber line at 
the leading edge; clockwise is 
positive. 


THETT is the camber angle at 
the trailing edge in degrees. 
It is measured from the chord 
line to the mean camber line at 
the trailing edge. Clockwise 
is positive. 


For blades whose chordline lies 
outside the blade profile, such 
as the DCA blading discussed in 
this report, extra grid lines 
surrounding the blade are 
needed to interpolate the blade 
position. The truncated value 
of MS( )/CAMPER is added to 

MS( ). The maximum allowed 
value of MS( ) + MS( )/CAMPER 
is 30 grid lines. These are 
the grid lines that enclose the 
blade profile. 


STABAC is used only for test 
cases where no blade shape is 
to be input. Default value is 
usually adequate. 
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CHOP Real 0.99 If leading edge or trailing 
edge radius is less than 2% of 
chord consult Ref. 9 for 
clarification. Normally, 

0.9 < CHOP < 1.0. 


SMOOTH Logical False For automatic addition of more 
blade definition points in the 
region of the leading and trail- 
ing edges set SMOOTE = .TRUE. 
This should be done for all 
blades except cusps and wedges. 


LEONLY Logical False LEONLY = .TRUE. for smoothing 
about the leading edge only; 
this is used if the trailing 
edge is a cusp or a wedge. 
SMOOTH must also be .TRUE. 


SLP1 Real 1 These parameters control the 

SLP2 Real 2 concentration of grid lines, if 

SLP3 Real 1 desired. The default values 

SLP4 Real 1 worked well for the DCA blades 
reported herein. For con- 
trolling the amount and loca- 
tion of the concentration 

we consult Ref. 9. 


The following logical and case control parameters are 


required only for flow solution runs. 


Name Type Default Description 
MINF Real None This is the Mach number at the 


upstream boundary. This can be 
determined from the non- 
dimensional velocity X output 
from the cascade tunnel data 
reduction program ‘Redd 5' and 
the relationship, 


Ys 1 y2 = x? which yields 
2 1 - x? : 


2 45 


Sires re. 
1-xe- Y-1 
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TOLS(I) is the tolerance for 
MS (I) and NOZES(I). There are 
three forms of input permitted. 


A) -1.0 < TOLS(I) < 0.0: Cal- 
culations of the flow solu- 
tion will proceed until the 
relative circulation error 


N 
© carc ~ Cexacr 
Cexact 


< |TOLS (I) | 


TOLS values between -107? 
and -10°° are typical 
values for grids. This 
method of input is appro- 
priate only for lifting 
(non-symmetric blades) 
cases. 


B) 0.0 < TOLS(I) < 1.0: Cal- 
culations of the flow solu-~ 
tion will proceed until the 
average relative change in 
potential is less than the 
absolute value of TOLS(I). 


(5) < |ToLs (I) | 
% /aVE 


Typical values should be 


3 and 107°. 


between 10° 

Cc) 1.0 < TOLS(I): Calcula- 
tions proceed until the 
number of iterations equals 
TOLS (I). 


Regardless of TOLS input, the 
solution for each grid will 
stop after 300 iterations if 
the TOLS criteria has not yet 
been made. All three forms 
were used for the cascade con- 
figuration reported herein 
with no discernible differences 
in results. 
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OVEREL Real 
UNDERL Real 
SUPREL Real 
NOWREL Integer 
NOTYET Integer 
TEGARD Real 
DAMP Real 
CII Real 
It Integer 
ALLOUT Logical 
we 
QUASI3 Legical 
NSTRM Integer 


\ 
| 


SCrNNNEFE eH 
Oe ee 


ee 
NOOO 


10 


- FALSE. 


- FALSE. 


Glas i aaa, alee aaa tite ne eee ne oe Ss een a ee 


The default values of these 
parameters are adequate for 
flow conditions in the subsonic 
cascade wind tunnel. 


Number of iterations between 
intermediate printouts of re- 
Siduals and Mach number. The 
information controlled by this 
parameter is of limited value 
in comparing with measured 
data, so a value greater than 
10 reduces the amount of com- 
puter printout. For the study 
reported herein 40 was used. 


To list the flow quantities at 
all grid points in the last 
mesh set ALLOUT = .TRUE.. Un- 
less a very coarse grid is 
used, the output resulting from 
ALLOUT = .TRUE. would be ex- 
tremely voluminous and of 
limited value. Until the cas- 
cade is configured so it is 
possible to take data from 
between the blades, ALLOUT 
should be .FALSE.. 


QUASI3 = .TRUE. to activate 
streamchannel thickness and/or 
radius variations. The cascade 
wind tunnel has no radius 
variations, but to simulate 3-D 
effects the streamchannel 
thickness is reduced at the 
exit boundary by a factor of 
1/AVDR. This data is placed in 
a file used by QSONIC if a 
quasi 3-D solution is desired. 


This is the position of desired 
streamsurface data on the 
streamchannel file used if 
QUASI3 = .TRUE.. Currently the 
default value of 1 identifies 
the proper streamsurface data 
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RINF Real 


ie 


WAKE Real 
MINF2 Real 
OMEGA Real 
VAXIAL Real 
FLOCO Real 


1 


0 


in the streamchannel file. If 
the output from the NASA code 
'Meridl' is used for the 
streamchannel data, then by 
using different values of 
NSTRM, different streamsurface 
data may be used. 


This is the spanwise radius 
at the upstream boundary di- 
vided by aerodynamic chord. 
Radius effects are activated 
if RINF # 1.0. The current 
version of QSONIC allows the 
following cases. 


QUASI RINF Results 
1) .FALSE. 1.0 Planar 2D flow 


2) .TRUE. #1.0 Thickness on 
file; radius 
on file. 


3) .TRUE. 1.0 Thickness on 
file; con- 
stant radius. 


Only 1) and 3) apply to the 
cascade wind tunnel. 


These parameters apply only if 
the test cascade is rotating 
and/or the downstream Mach is 
near !.0. 


At this point all of the logical and case control para- 


meters necessary to use QSONIC for the flow conditions pos- 


sible in the subsonic cascade wind tunnel have been discussed. 


The following is a description of the Bulk Data input for 


the geometry generation run and the flow solution run. 
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The format for all variables in the bulk data for mesh 
generation (geometry) is namelist form. The namelist is 


INSTUFF. 


Name Type Default Description 


H2 Complex None This is a table of points de- 
fining the blade profile. The 
real vart = X, and imaginary 
part = Y coordinate. The table 
begins at the point of maximum X 
value at the trailing «dge and 
proceeds clockwise bacx around 
to the first point, which is re- 
peatec. The blade must be at the 
Stagger angle and the origin at 
the point of minimum X. For 
Electrostatic Analog grids, the 
Stagger angle must be positive 
(leading edge low, trailine edce 
high). The maximum number of 
points in H2 is 99 for the inter- 
polation scheme or 63 for the 

we Electrostatic Analog. 


BUG2 Logical -FALSE. BUG2 = .TRUE. for a more detailed 
output of geometry generation. 
This will include the X,Y coordi- 
nates that define the mesh as 
well as second derivatives at 
grid points on the body. Except 
for trouble shooting this data is 
of limited value at the vresent 
time since there is no graphic 
Output. 


The bulk data for flow solutions consists of a mesh 
file and the streamchannel data file if quasi-3D effects are 
to be calculated. The mesh file is created by QSOMNIC during 


the mesh (grid) generation run. No further inputs are re- 


quired from the user for the mesh file. 


144 


The streamchannel data file must contain a table of 
streamtube thicknesses, radial positions and corresponding 


X values along the streamsurface. 


Name Type Default Description 
CHO Real None CHO is the aerodynamic chord 


multiplied by the cosine of 
LAMDA0. (LAMDAO = stagger angle) 


NRSP Integer None NRSP is the total number of data 
points in each of the tables of 
thickness, radial position, and X 
location. If MRSP is 2, a linear 
distribution is obtained between 
the endpoints given. NRSP = 2 
was used for the study reported 
herein with good results. 


RM Real None Array of corresponding X loca- 
tions for thickness and radius 
data values. X = 0 represents 

the leading edge of blade, with 
we the blade at stagger angle. The 
units can be any consistent 
length scale common to.CHO, RM, 
RMSP and BESP. Inches were used 
in this study. 


RMSP Real None Spanwise radial positions of 
streamsurface at the X locations 
given in RM. RMSP was not used 
in the current study. 


BESP Real None Array of streamtube thickness 
values at the X locations speci- 
fied in RM. For the study re- 
ported herein, at X=0Oa 
streamtube thickness of 1.0 was 
arbitrarily selected. The 
streamtube contraction through 
the test section was similated 
by reducing the thickness at 
X = 0, by a factor of 1/AVDR at 
the trailing edge. Duval [Ref. 
3] explains AVDR. 
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D.3 PREPARING INPUT FILES 
QSONIC was originally configured to use several input/ 
output devices while reading data and generating output. 


The input/output devices are listed below as used by QSONIC, 


I/O Unit ‘Usage 
2 File containing streamchannel data. This is used 
only if QUASI = .TRUE.. 
5 Standard card input; NAMELISTS PARAMS, INSTUF. 
6 Standard printed output. 
13 For mesh generation runs, coordinates of all mesh 


points are written here. For flow solutions, xX, 
Y, velocities, pressures and MACE are recorded 

for graphic display. The program currently has no 
graphic output capability. No user action is 
necessary to create this file. 


18 Used as temporary storage. No output is stored 
here. No user action is necessary in conjunction 
with this file. 


23 Previously developed mesh coordinates are read in 
from this file during the flow solution. After a 
mesh generation run, the user must create this 
file and put in it the data from I/O unit 13, so 
that during the flow solution QSONIC can read in 
the mesh points. 


QSONIC is presently configured to operate with the CMS 
system of the IBM 370 computer. This system provides a high 
degree of flexibility in parameter selection. With this 
‘es system, all the input/output units previously mentioned are 
on the disk space assigned to the Turbopropulsion Laboratory 
(TPL). Access to QSONIC and the TPL disk space can be ob- 


Le tained through the Director of TPL. 
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The first step in using QSONIC is the creation of the 
data file necessary for the mesh generation run. ‘This is 
done using the standard procedures of the XEDIT function of 
the CMS operating system. Reference 14 has specific in- 
structions for creating new files. The filename and file- 
type for the data used in this study was NAMELIST™T GFOMD. 
Once the data file is opened, the necessary data is input 
beginning at column 2 of the virtual card. Since the varia- 
bles are in namelist form, they can be input in any order. 

Table D.1 is an example of the data file necessary for 
the mesh generation run. The TITLE must appear on the first 
line (FORMAT = 20A4). After the TITLE, the logical and case 
control parameters are input after the namelist GPARAMS. 
When all the case control parameters required are input the 
PARAMS namelist is closed with &END. On the next line of 
the data file the buik data for the mesh generation run is 
input in namelist form, with the namelist &INSTUF. The 
Hl = 100*(0.0,0.0) that appears after &INSTUF on Table D.1 
was used on earlier versions of QSONIC, but is not used in 
the present version. It should, however, appear in the data 
file before the H2 variables (X,Y coordinates defining the 
blade profile). 

At this point, some discussion of the coordinates de- 
fining the blade profile is warranted. Table D.2 is a 
listing of the X and Y coordinates of the DCA blading used 


in this study. Figure D.1 is a plot of these coordinates. 
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Recall that the coordinates defining the blade profile for 
QSONIC must be for the blade at the stagger angle and mini- 
mum X at the origin. The coordinates of Table D.2 were 
translated ane rotated using a coordinate transformation 
routine for the HP-67 programmable calculator. These new 
coordinates appear in the namelist INSTUFF on Table D.1. 
Figure D.2 is a plot of the translated and rotated coordi- 
nates. It is highly recommended that such a plot be made 
for any new blade profiles, to ensure that the original 
coordinates are translated and rotated properly. 

The second step in using QSONIC is the creation of the 
Gata files necessary for the flow solution run. The file 
used for the flow solution in this report is on the TPL disk 
space with a filename/filetype of NAMELIST FLOYD. The sim- 
plest way to open this data file is to use the YEDIT func- 
tion, as discussed in Ref. 14, to start a new file. Then 
input the same data as is in the data file for the mesh 
generation run usine the X™DIT subcommand GFT (filename) 
(filetype). The appropriate changes anc additions can then 
be made to this file. Table D.3 is an example of the data 
file just discussed. 

Two more data files are recuired for the flow solution. 
The data for one of these is created by the mesh generation 
run. ‘The other file contains the streamchannel thickness 


data for implementing quasi-3D effects. 
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After the mesh generation run, a file with the file- 
name/filetype MESHGEN DATA will appear on the disk. Create 
a new file with the filename/filetvne MESHIN DATA. This is 
most easily done by issuing the command 'XEDIT MESHIN DATA'; 
then use 'GET MESHGEN DATA'. This file contains the pre- 
viously developed mesh coordinates. 

The streamchannel data file should have the filename/ 
filetype of DATA5SD DATA. The format for the data file is 


shown below. 


Virtual Card Column No. Variable Name 

L BLANK 

2 BLANK 

3 21-30 CHO 
4 BLANK 

5 36-40 NRSP 
6 BLANK _- 

7 BLANK 

8 

9 
10 
11 
1-80 (8F10.5) RM 


As needed 


1-80 (8F10.5) RMSP (not used in 
this study) 
As needed 


1-80 (8F10.5) BESP 


Table D.4 is an example of the data file used in the 


present study. Since NRSP = 2 was used for this study a 
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linear distribution is assumed for the streamtube thickness 
values and only 2 values of RM and BESP are required; there- 


fore, only 1 virtual card was required for each array. 


D.4 PROGRAM OUTPUT 

The output generated by QSONIC for the geometry genera- 
tion run includes a printed listing (I/O unit 6) and a mesh 
point file (I/O unit 13). The printed listing under the CMS 
system I/O unit 6 is normally the computer terminal unless 
the command 'FILEDEF 06 PRINTER' has been invoked. It is 
unusual for the program to run properly the first time, so 
initially it is helpful to have the printed listing appear 
at the terminal. Once the program is running properly the 
output should be sent to the line printer. 

The flow solution run output consists of a printed 
listing (I/O unit 6) and a plot data save file (I/O unit 
13). Once the flow solution is running properly the printed 
listing should be sent to the line printer. 

Table D.5 is an example of the output generated by the 
geometry generation run. Figure D.3 is a plot of the grid 
output points on the blade surface, horizontal chord, pro- 
duced by the mesh generation run. Figure D.4 is a plot of 
the grid output points with the blade at the stagger angle. 
Tabie D.6 is an example of the output generated by the flow 


solution. 
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A det. ied explanation of the printed output for the 


program QSONIC may be found in Ref. 9. 


D.5 RUNNING THE PROGRAM 
The files on the TPL disk space that apply to OSONIC 

are listed below: 

QSONIC EXEC Al 

QSONIC FORTRAN Al 

QSONIC TEXT Al 

NAMELIST GEOM Al 

NAMELIST FLOW Al 

NAMELIST GEOMD Al 

NAMELIST FLOWD Al 

DATAS DATA Al 

DATA5D DATA Al 

QSONIC EXEC sets the input/output devices required to 

read and store data. QSONIC FORTRAN is the source program. 
To document the changes to QSONIC necessary to use the code 
with the IBM 370 operating system and serve as a reference 
for future users, a program listing is included at the end 
of this appendix. QSONIC TEXT is the computer executable 
code created when QSONIC FORTRAN is compiled. NAMELIST GEOM 
and NAMELIST FLOW are the data files for the geometry gen- 
eration and flow generation respectively for the example in 
Ref. 9. DATAS DATA is the streamchannel data required for 


the quasi-3D solution for the example in Ref. 9. 
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NAMELIST GEOMD is the data file for the ceometry gen- 
eration for the DCA blading used in the studv reported 
herein. NAMELIST FLOWD is the file for the flow solution 
for the cascade configuration used in this study. 

QSONIC expects the input data to be in a file on the 


TPL disk space named NAMELIST DATA. Since the first time 


QSONIC is run is to develop the body centered mesh, the file 
NAMELIST GEOMD must be renamed NAMELIST DATA, using vroce- 
dures specified in Ref. 13. Because QSOMIC reauires large 
amounts of virtual memory, extra storage must be defined for 
the code to operate. This is accomplished by issuing the 
command 'DEFINE STORAGE 1504K'. 

With the data file renamed and more storage defined, 
type 'QSONIC' to load the program. The output will appear 
on the terminal screen unless FILEDEF $6 PRINTER was invoked 
prior to loading the program. 

After the mesh generation is complete, rename NAMELIST 
DATA to NAMELIST GEOMD and change NAMELIST FLOWD to NAMELIST 
DATA. Create a data file with the filename/filetype MESHIN 
DATA. The elements of this file are the same as the elements 
in the file MESHGEN DATA that was created by the mesh gen- 
eration run. The necessary input/output files are now con- 
figured for a flow solution run. Issue the command 'OQSONIC’ 
to begin execution. 

If the program output appears at the terminal it is 


possible to have some I/O error messages appear with the 


output. This is because the write statements in QSONIC are 
formatted for the 132 character long line of the printer. 
These errors do not affect the validity of the program 
output. 

The explanation for any error or condition message 


generated by QSONIC can be found in Ref. 9. 


D.6 QSONIC UPDATE 

Recently an improved version of QSONIC was reported by 
NASA Lewis Research Center [Ref. 13]. The new version re- 
quires less virtual memory and executes approximately 303% 
faster than the version presently in use at NPS. Also, the 


output appears in a different format than is described in 


this appendix. Reference 18 describes the most recent ver-— 


sion of QSONIC in detail. 
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INPUT DATA FOR MESH GENERATION RUN OF QSONIC 
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